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The  significant  role  played  by  microdefects  (damage)  in  the  process  of 
crack  formation  and  growth  is  commonly  recognized.  We  distinguish  two  extreme 


cases  of  the  influence  of  defects  on  the  fracture  process;  modeling  of  these 
two  cases  requires  essentially  different  formalisms. 

Case  1  -  The  intensity  of  damage  formed  as  a  response  to  the  stress  concentra¬ 
tion  at  the  tip  of  a  propagating  crack  is  much  greater  than  the  intensity  of 
the  pre-existing  damage.  The  crack  propagation  is  then  Inseparable  from  the 
evolution  of  the  damage  accompanying  the  crack.  This  strongly  cooperative 
phenomenon  is  modeled,  for  example,  by  Crack  Layer  Tbeoryl  based  on  thermo¬ 
dynamics  of  irreversible  processes. 

Case  2  -  A  crack  propagates  through  a  pre-existing  field  of  defects  causing 
negligible  changes  to  the  field.  Hie  fluctuation  of  the  microdefect  field  is 
directly  reflected  in  stochastic  features  of  fracture  surfaces  and  also  leads 
to  the  scatter  of  experimentally  observed  fracture  parameters,  such  as  criti¬ 
cal  crack  length,  critical  load,  etc.  A  probabilistic  approach  seems  most 
adequate  under  these  circumstances. 

Below,  we  present  a  probabilistic  model  of  crack  formation.  Our  consid¬ 
eration  is  based  on  the  approach  outlined  earlier2. 

II.  Probabilistic  Model  of  Crack  Formation 

We  address  the  following  problem:  what  is  the  probability  of  crack  for¬ 
mation  beyond  a  particular  depth  in  an  elastic  medium  with  fluctuating 
strength  properties  under  prescribed  loading  conditions. 

Let  us  consider  an  ensemble  or  macroscopically  identical  specimens  sub¬ 
jected  to  Identical  loading  conditions.  Let  the  resulting  fracture  profiles 


be  superimposed  on  the  same  plot  as  shown  In  Fig.  1  (such  experiments  are 
reported3 »** 5  f0r  simple  tension,  and®  for  four  point  bending).  Apparently, 

each  crack  trajectory  is  unique  (no  two  coincide).  More 
importantly,  in  a  given  test  from  the  above  sequence  of  identical  tests,  any 
of  the  plotted  trajector  for  each  specimen  ies  must  be  viewed  as  a  priori  pos¬ 
sible.  Statistical  analysis  of  the  observed  trajectories  allows  one  to  formu¬ 
late  reasonable  assunptions  about  the  nature  of  the  set  0  of  all  possible 
crack  trajectories  for  each  specimen  under  given  test  conditions.  Thus  Fig.  1 
presents  a  sample  subset  of  0. 

If  we  assume  that  only  one  crack  is  formed  in  eacH  specimen  (single  path 
fracture),  then  the  probability  P(X)  that  a  crack  is  formed  beginning  at  the 
notch  and  extending  to  or  beyond  the  depth  X  can  be  written  as 

P(X)  =£  P{X|u}P{u)}  .  (1) 
0 

where  P{w}  is  the  probability  that  the  crack  'chooses'  a  path  w  among  all 
possible  paths,  and  P {X I  id)  is  the  conditional  probability  that  the  depth  X 
is  reached  by  a  crack  formed  along  w.  In  a  continuum-based  model,  the  space 
fl  is  uncountable,  so  the  sum  in  (1)  has  to  be  substituted  by  an  integral 

P(X)  -fra  Ioj)  d  (w)  .  (2) 
Q 

This  brings  about  three  tasks:  to  propose  an  adequate  crack  path  space  0, 
to  choose  a  probability  measure  dp(w)  on  0,  and  to  determine  the  probabil¬ 
ity  P  {X I  id)  . 

The  first  two  tasks  are  closely  related.  Their  solution,  in  principle, 
is  based  on  a  statistical  analysis  of  fracture  surfaces,  and  preliminary  data 
has  been  obtained  recently.  Here  we  consider  a  'diffusion'  model,  i.e.. 


we 


3 


model  crack  trajectories  by  Brownian  paths  y  =  u(x),  0  _<  x  <  X,  w(0)  =  0,  (t^ 

7 

*=  {o)J  will  denote  the  space)  and  choose  dp(«)  to  be  a  Wiener  measure 


(w)  =  const  •  exp 


[a) '(x) ]2  dxj  fl  dw(x) 
Jx=0 


(3) 


Below  we  will  meet  the  parameter  D  >  0  from  (3)  in  a  diffusion-type  equation 
(13)  for  a  probability  of  crack  formation.  By  analogy,  we  refer  to  D  as 
'crack  diffusion  coefficient'.  D  reflects  the  tendency  of  crack  trajectories 
to  deviate  from  the  X-axis  and  is  experimentally  measurable5.  Evaluation  of 
D  uniquely  determines,  in  the  'diffusion'  model,  the  set  of  all  possible 
crack  trajectories  as  well,  as  the  measure  in  (2). 

The  third  task,  determination  of  the  conditional  probability  PtXlu),  is 
based  on  the  formulation  of  a  fracture  criterion. 

III.  Probability  of  Crack  Formation  Along  Given  Path 

In  the  case  of  brittle  fracture  of  a  homogeneous  material .  the  Griffith 
criterion  of  uncontrolled  crack  propagation  is  commonly  written  as 


IOJ  I 

=  2r  .  —  >  0  '  (4) 

z  =  z  '  z  =  z 

c  c 

where  J  is  the  energy  release  rate,  r  is  the  specific  surface  energy,  and 

Zc  is  the  critical  crack  length  for  the  crack-loading  configuration  under  con¬ 
sideration.  Strictly  speaking,  (4)  is  a  condition  of  global  instability,  if 

dJ/dZ  >  0  for  all  l  >  Namely,  (4)  then  assures  that,  as  the  crack 

advances,  there  is  sufficient  amount  of  potential  energy  released  to  compen¬ 
sate  for  the  energy  required  for  new  surface  formation 


4 


J<£)  >  2r 


(5) 


for  2.  >  2C 

In  a  heterogeneous  material,  y  becomes  a  random  field  which  we  assume 
to  be  statistically  homogeneous.  Now  we  take  (5)  as  the  criterion  of  local 
crack  propagation.  Thus,  for  a  heterogeneous  strength  field  y,  the  global 
instability  criterion  should  be  written  as  an  integral  requirement  that  (5)  is 
met  everywhere  along  a  fracture  path.  We  further  assume  that,  at  distances 
exceeding  certain  r,  values  of  the  field  y  are  independent.  This  correla¬ 
tion  distance  r  will  be  assumed  to  be  much  smaller  than  the  crack  size. 

Subdividing  the  crack  trajectory  into  small  portions  and  assuming  that 
the  crack  arrest  at  each  step  is  an  unlikely  event  with  the  probability  pro¬ 
portional  to  the  length  of  the  step,  we  arrive  at  the  following  expression  for 
the  probability  of  crack  penetration  to  a  depth  X  along  a  particular  path  <d 
(see  Appendix) : 


P(XM 


exp 


P{2r  >  JJx)) 


(6) 


where  P(2y  >  Jw(x))  stands  for  the  probability  that  the  value  of  the  random 
field  y  exceeds  the  value  of  the  energy  release  rate  at  a  current  point  (x, 
w ( x) )  of  the  path  u. 

To  describe  P  (2y  2.  <Jw(x)},  we  assume  that  locally  a  crack  chooses  the 
'easiest'  path,  i.e.,  it  minimizes  the  energy  barrier  2y  -  J.  This  is  equiv¬ 
alent  to  the  assumption  that  y  along  crack  trajectories  represents  minimal 
values  of  the  y-field,  since  fluctuations  of  the  y- field  occur  at  a  scale 
small  in  comparison  to  crack  increments,  considered  for  energy  release  rate 

g 

evaluation.  According  to  the  theory  of  statistics  of  extremes,  we  employ  a 
Weibull  distribution  for  y  at  every  point  of  crack  trajectory 
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F(t> 


1  -  exp 

0 


(^r)‘ 


(*  2  *min) 
(*  <  ^min* 


(7) 


where  a  >  0,  rmin  and  (commonly  kncwn  as  shape,  minimal  value  and  scale 

parameters)  are  empirical  constants.  Then 


p{2T  >  J„(x))  =  1  -  F(Ju(x)/2)  =  exp 


Ymln^ 


(8) 


Putting  Eqs.  (2),  (6)  and  (8)  together  we  get  the  following  expression 
for  the  probability  that  a  crack  is  formed  to  a  depth  greater  than  X: 


(9) 


where  d^  (w)  is  given  by  (3). 

7 

The  functional  integral  in  (9)  can  be  evaluated  by  various  techniques  , 
some  of  which  we  use  below. 


IV.  Crack  Diffusion  Equation 

In  this  section  we  reduce  integration  in  (9)  to  solving  a  conventional  dif- 

9 

fusion  type  equation.  This  is  achieved  by  representing  P(X)  as 

Oo 


P(X) 


P(X, Y)dY 


(10) 


—  &G 


where  the  new  function  P(X,Y)  is  introduced  by 


P{2r  >  V*)) 


dx 

r 


dM<»>(u» 


(11) 
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$ 


Here  P{2y  >  J  (x)}  is  given  by  (8),  is  the  set  of  all  possible  crack 

trajectories  beginning  at  the  origin  and  crossing  a  point  (X.Y),  and  dpADy(«i>) 
is  the  corresponding  conditional  Wiener  measure.  The  meaning  of  P(X,Y)dY  is 
the  probability  that  a  crack  is  formed  beginning  at  the  notch  and  passing 
through  the  'window'  Y  .<  y  _<  Y  +  dY  at  the  depth  X. 

The  probability  P(2y  >  J^x)}  given  by  (8)  is  a  functional  of  *t.  since 
the  energy  release  rate  J^x)  is  a  functional  of  the  crack  trajectory.  A 
special  feature  of  this  functional  is  that  it  is  insensitive  to  small  pertur¬ 
bations  away  from  the  crack  tip,  l.e.,  functional  derivatives  such  as 

^Ja>(x)/6w  can  be  approximated  as  being  proportional  to  {-functions  coneen- 
triated  at  the  crack  tip.  It  follows  by  means  of  a  functional  Taylor  expan¬ 
sion  near  a  conditional  average 


u>0(x)  =  f  w(x)  dMx?Y(b,) 
°X,Y 


that  P{2y  >  <Jw(x)}  can  be  substituted  by  a  function  V  of  current  crack  tip 
coordinates  x  and  w(x)  thus  allowing  to  rewrite  (11)  as 

!j.DJ(u>)  (12) 

(  •'o  ) 

°X,Y 

Exploiting  results  of  M.  Kac®  concerning  integrals  of  the  form  (12),  we  con¬ 
clude  that  P(X,Y)  is  the  solution  of 


i 

P(X ,  Y)  = 

/  expj- 

f  V(x,(i>(x)  )dx  \ 

« 

J  (  ' 

J  0  * 
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P(X ,  Y) — ►  0  ,  as  T  —►±00  (X  >  0) 

P(O.Y)  =  5 (Y) 


(14) 


We  refer  to  (13)  as  'crack  diffusion  equation'. 

We  remind  that  P(X,Y)  gives  rise  to  the  function  P(X)  via  (10),  and 
consequently  determines, the  probability  density  p(X)  of  a  crack  penetrating 
e.:  JLy  to  a  particular  depth  X 


p(X) 


dP(X) 

dX 


(15) 


(by  its  definition,  l-P(X)  is  the  distribution  function  for  the  X-coordinate 
of  the  crack  tip).  Apparently,  P(X)  contains  less  information  than  P(X,Y). 
However,  it  allows  one  to  reconstruct  the  y-field  parameters  on  the  basis  of 
experimental  data,  as  discussed  in  section  VI. 


V.  Illustrative  Example 

In  this  section,  we  consider  the  effect  of  the  parameters  and  of 
the  crack  diffusion  coefficient  on  the  crack  depth  distribution  p(X)  for  a 
particular  example  of  specimen-loading  configuration. 

Let  us  consider  a  random  crack  formed  in  a  semi-infinite  plane,  under  the 
action  of  a  pair  of  forces,  as  shown  in  Fig.  2. 

In  Section  III  the  evaluation  of  the  probability  of  crack  penetration  to 
a  depth  greater  than  X  was  reduced  to  the  evaluation  of  the  energy  release 
rate  Jw(x)  as  a  functional  of  u)  (see  (9)).  The  energy  release  rate  dependency 
on  a  crack  trajectory  is  quite  complex;  no  general  solution  is  available. 
There  are  'long  wave'  and  'shortwave'  u- contributions  to  J^(x).  The  'long 
wave'  effect  is  accounted  for  by  perturbation  methods  or  other 
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means.*®' 11 ' The  'short  wave'  w-contribution  can  be  approximated  by 
existing  kinked  crack  solutions,11 '**  since  the  functional  is  highly  sen¬ 
sitive  to  the  direction  of  the  crack  at  the  very  tip  and  is  practically  not 
affected  by  small  perturbations  away  from  the  tip.  It  follows  from  an  analy¬ 
sis  of  the  existing  solutions  that  the  'short  wave’  effect  dominates.  Consid¬ 
ering  the  nature  of  the  kinks  in  the  'crack  diffusion  model',  we  take  the  fol¬ 
lowing  approximation  for  the  energy  release  rate  associated  with  the  crack 
formed  along  «•>  to  a  depth  x  : 

J^fx)  =  JQ(x)  -  k  JQ(x)  (<i )/r)2  (16) 

where  JQ(x)  stands  for  the  conventional  energy  release  rate  for  a  rectili¬ 
near  crack  of  the  length  x  formed  along  the  X-axis 


J0(x) 


2.2  71 

E  •  x 


(17) 


the  factor  k  =  0.5  is  evaluated  on  the  basis  of  a  nimerical  solution  for  a 
kinked  crack. 13 

Before  proceeding  to  the  evaluation  of  the  probability  P(X)  (see  (9)), 
it  is  convenient  to  introduce  the  mean  value  y*  of  the  y-field  instead  of 

Y0  in  the  Weibull  distribution  (7),  obtaining 


rv, 


1  -  exp 


T'min 

^min 


(T  >  Tmin* 


F(y)  = 


(18) 


(y  <  Ymin) 


We  used  (a  straightforward  calculation) 


oo 


y*  =  *min  +  c<o)*0  ' 


c(a)  =  /  exp 

0 


a 


j  dz 


(19) 
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The  corresponding  changes  in  (8)  and  (9)  are  immediate. 

Zeroth  Approximation  (indicated  below  by  the  subscript  0).  Substituting 
J^x)  in  (9)  by  Jc(x)  (see  (16,17)),  and  performing  the  functional  inte¬ 
gration  (which  is  trivial  since  the  integrand  does  not  depend  on 

a)  )  ,  we  obtain 


P0(X)  =  exp 


/'  X  r 

{-/  ' 

'  0  1- 


c(a) 


J0(x)/2  -  Tmln' 


r*  ~  Tmin 


“) 


(20) 


It  is  convenient  to  rewrite  (20)  in  a  dimensionless  form. 


Let  us  introduce  a  'Griffith  length'  £Q  determined  from  JQ(^G)  =  2y* 
and  representing  the  length  of  an  equilibria!  rectilinear  crack  along  the 


X-axis  in  an  ideally  homogeneous  material  with  the  fracture  energy  y*.  This 
yields  dimensionless  crack  penetration  depth  and  correlation  distance 

z  =  x/£g  ,  p  =  r/£Q  .  (21) 

We  also  have  from  (17) 


1 

-  =  —  =  -  .  (22) 

2y*  x  z 

Passing  to  the  dimensionless  quantities  in  (20)  ,  we  get 


where 


q  =  W r*  (0  <  q  <  1) 


(23) 


(24) 


The  probability  density  pQ( Z)  =  -  dPQ(Z)/dZ  of  the  dimensionless  crack 
penetration  depth  Z  fcf.  (15))  is  shown  in  Fig.  3  for  various  q  and  a  = 


3 


i-2  15 


10 


Notice  first  that  the  mean  crack  depth  <  X  >Q  «  Jig  <  Z  >Q  does  not 

coincide  with  the  depth  Iq  of  the  'Griffith  crack*  corresponding  to  the  mean 
value  of  the  fracture  energy  (as  should  be  expected). 

Secondly,  as  q  increases,  the  mean  crack  depth  increases,  the  variance 
of  the  crack  depth  decreases  (though  non-mo notonically  in  general),  and,  when 
q  approaches  1,  P0(Z)  tends  to  S(Z-l).  The  latter  is  a  consequenoe  of  the 

fact  that  Ymin/y*  =  1  corresponds  to  the  distribution  F(y)  being  a  step- 
function,  i.e.,  the  field  y  being  a  homogeneous  non-random  field  with  the 
value  y*. 

Our  final  remark  on  the  zeroth  approximation  concerns  a  straightforward 
way  to  arrive  at  (20)  without  functional  integrals.  Namely,  let  us  consider  a 
rectilinear  crack  propagating  along  the  X-axis  through  the  same  random  field 
y,  as  above.  Then  adopting  the  Poisson  process  type  assumption  (see  Appen¬ 
dix,  specially,  (A7))  one  arrives,  in  a  standard  way,  at  the  differential 
equation 

/  dPc(X)  j 

—  =  ;P{2^  V*)}  *  p0w 


(25) 


P0(I) 


X  =  0 


=  1 


whose  solution  is  exactly  (20)  (in  view  of  (8),  (17)). 


Next  Approximation.  Now  we  consider  Jw(x)  given  by  (16).  As  in  the  zeroth 
approximation,  it  is  convenient  here  to  pass  to  the  dimensionless  quantities 

z  =  x/£G<  z  =  x/Zq,  p  =  r/  g,  q  =  and*  in  action,  <?(z)  =  u>(z)/r. 
Substituting  (16)  into  (9),  using  (22).  and  taking  into  consideration  a  simple 
relationship  for  functional  integrals 
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/  (...) 


dMxD)U> 


/(...) 


0 


X 


„  JA),  .  4  •  D 

,)  d/xz  (*)  a  =  _ 

°z  '  r2 


(26) 


we  get  the  following  expression  for  the  probability  of  crack  penetration  to  a 
depth  greater  than  Z: 


P(Z) 


/eipl:-C 


exp 


-  o(a) 


(l-k^>2)/z-q\  1  dz  ) 

1  -  q  ;  J  p  } 


d/z>(^)  (27) 


Maintaining  in  the  inner  integrand  the  terms  of  the  degree  up  to  two  (in  <e  ) , 
we  obtain 


P(Z)  =  PQ(Z)  •  £  exp  (■/:  u(z)*2(z)dz j  (<p)  , 


(28) 


where  PC(Z)  is  the  probability  of  crack  pentration  to  a  depth  grater  than  Z 


evaluated  in  the  zeroth  approximation  (23); 

ok  / c(a)  \  (i/z_q) 

u(z)  =  - l— —  I  -  exp 

P  \  1-q/  z 


-  c(a) 


1/z-q 

1-q 


(29) 


Evaluation  of  the  functional  integral  in  (33)  can  be  reduced  to  solving 
an  ordinary  differential  equation,  1*  namely. 


J "  exp  j-,/*  u(z)v>2(z)dz  j  d^z^  (<?)  =  [Q(z)l 


(30) 


]z  -  0 


where  the  function  Q(z)  satsfies 


Q  (z)  -  A’u(z)  -Q(z)  =0 

Q(z)|  =1  , 

'z  =  Z 


(31) 


Q'(z)  =  o 

z  *=  Z 


The  probability  density  p(Z)  =  -dP(Z)/dZ  of  the  dimensionless  crack  penetra¬ 
tion  depth  Z  (cf.  (15))  is  shown  in  Fig.  4  for  various..,  q  =  0  and  a,  p  as  above. 
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Notice  that,  as  D  increases,  the  mean  crack  depth  decreases.  This  can  be 
understood  in  the  following  way.  In  terms  of  energy,  one  finds  lower  average 
energy  release  rate  for  trajectories  associated  with  higher  D,  since  those 
have  steeper  kinks.  In  addition,  the  statistical  weight  of  such  trajectories 
increases  with  D. 

Due  to  this  effect,  evaluation  of  y*  on  the  basis  of  measured  average 

2 

crack  depths  (y*  =  1.1?"  /(E-<  X  >))  ignoring  'crack  diffusion'  would  result 
in  overestimation  of  the  strength.  More  generally,  reconstruction  of  the 
y-field  parameters  requires  crack  diffusion  analysis.  The  latter  can  be  based 
on  a  straightforward  statistical  analysis  of  fracture  surfaces  in  an  ensemble 
of  identical  specimens5  or  on  a  more  sophisticated  analysis  of  a  single  frac¬ 
ture  surface  from  a  sufficiently  large  specimen. 


VI.  Conclusions 

In  brittle  solids,  specific  fracture  energy  evaluated  on  the  basis  of 
fracture  experiments  displays  large  scatter  and  specimen  geometry  dependence. 
The  probabilistic  model  outlined  above  treats  specific  fracture  energy  y  as 
a  random  field  and  suggests  certain  experimental  procedure  for  evaluation  of 
the  y-field  parameters. 

Two  types  of  parameters  enter  the  expression  for  the  probability  distri¬ 
bution  of  crack  penetration  depth:  Weibull  distribution  parameters  (pointvise 
characteristics  of  the  y-field)  and  the  correlation  distance  together  with  the 
'crack  diffusion  coefficient',  both  reflecting  non-local  properties  related  to 
material  structure  such  as  grain  size.6  The  latter  ones  can  be  extracted  frem 
fractographic  analysis  of  a  sample  set  of  crack  trajectories.  Then  a  compari- 


son  of  the  mathematical  expectation,  variance  and  other  moments  of  the  crack 
penetration  depth,  as  rendered  by  the  model,  with  those  found  experimentally 
yields  the  values  of  the  Weibull  distribution  parameters. 

Thus  determined  y-field  parameters  can  be  used  (a)  to  rank  brittle 
materials  with  regard  to  toughness;  <b)  to  predict  structural  reliability. 

The  dimensionless  'crack  diffusion  coefficient*  A  =  ^ *  D/r^  introduced 
above  gives  a  new  similarity  criterion  for  single  path  fracture  (to  be  used  in 
addition  to  similitude  in  geometry  and  stresses) .  It  involves  a  mixture  of 
macroscoplcal  and  microscopical  parameters,  thus  constraining  the  choice  of 
materials  for  small  scale  modelling  of  brittle  fracture. 


APPENDIX  -  CALCULATION  OP  P{X|w) 


We  subdivide  the  interval  0  <  x  <  X  into  m  equal  portions  Ax^  whose 
length  Ax  is  small  enough  so  that  can  be  considered  constant  over 

each  Ax^.  At  the  sane  tine,  we  assume  that  Ax  is  much  greater  than  the 
correlation  distanoe  of  the  y-field,  Ax  >>  r.  Let  us  further  split  each  Axk 

9 

into  n  *  Ax/r  subintervals  6xj  ((k  -  l)n  <  j  <  kn)  of  the  length  r  (see 
Fig.  5)  and  consider  y  to  be  constant  over  each  6xj.  ^us  the  interval 
0  <  x  <  X  is  subdivided  into  m  •  n  subintervals  6xj. 

Crack  formation  between  the  notch  and  a  depth  Z  means  failure  over  all 


8x  's.  According  to  the  energy  criterion  of  local  failure,  the  failure  over 

J 

Axj  occurs,  if  the  released  potential  energy  Allj  is  greater  than  the 


energy  required  for  new  surface  formation 


inJ  >  •  6*j-18 


Thus 


mn 

P{Xlw)  =  P  {  fl  <  Allj  >  2rj  •  Axj)  (  (A1) 

J-l 


To  evaluate  the  probability,  we  consider  crack  formation  as  an  ordered 
sequence  of  local  failures.  Therefore  it  is  convenient  to  represent  P{Xlu} 
as  the  following  product  of  conditional  probabilities: 


P{X  I  w]  *  11  P  {  (  AHj  >  2Tj  •  6Xj)  |  fl  (  ATIj  >  2Ti  •  Ax^  }  .  (A2) 

J-l  1  i=l 

The  potential  energy  released  at  the  J-th  step  under  the  condition  that  all 
preceeding  elements  failed  can  be  evaluated  through  the  conventional,  in  frac¬ 
ture  mechanics,  energy  release  rate  Jw(x) 


A°J  *J«o(xJ)*6^ 


Thus  (A2)  becomes 


P{Xlu>) 


n  p  {ju(Xj)  >  2Tj} 

j»i 

m  kn 

-  n  n  p  (jw(xj)  >  2Tj} 

k=l  j=(k-l)n+l 


(A4) 


where  the  inner  product  Represents  the  (conditional)  probability  of  failure 
over  Axk.  According  to  the  choice  of  Ax,  Jw(x)  is  (approximately)  constant 
over  each  Axk,  i.e.,  for  a  given  k,  all  the  values  Jw(xj)  in  the  inner 

i 

product  in  (A4)  can  be  substituted  by  the  one  and  same  value  J  (x,  ).  Also 

w  kn 

Tj's  are  independent  random  variables,  which  have  the  same  (Weibull)  distri¬ 
bution  (7).  We  conclude  that  the  probabilities  in  (A4)  are  all  equal,  for  a 
given  k.  Passing  to  the  opposite  event  in  each  of  the  n  =  Ax/r  (equal)  fac- 

I 

tors  in  (A4) .  assuming  that  its  probability  is  small,  and  using 

(1  -  e)D  »  exp  (-ne)  (small  e)  ,  (A5) 

» 

I  we  rewrite  (A4)  ,  as  follows: 


P{xlw)  =  n  (  1  -  P  {2y  >  Jw<xkn») 


Ax/r 


k=l 


a  exp 


Ax 


-r  f  (»r  >  •V’W  — > 

k=l 


( A6) 


J  Substituting  the  above  summation  by  integration,  one  arrives  at  (6). 

It  is  worth  noting  that  one  can  arrive  at  (6)  by  a  standard  Poisson  pro¬ 
cess  type  argument,  if  one  postulates  that  the  probability  of  the  crack  arrest 

i 

over  a  segment  x  <  x 1  <  x  +  Ax  is  proportional  to  Ax 
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P  {  2y(x'  ,«(x ') )  >  J  (x*>  for 

1  _  u) 

=  P{2y  >  JU(X)} 


some  x 


<  x '  <  x  +  Ax  | 


Ax 

r 


(A7) 
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A  ensemble  of  macroacopically  identical  specimens  cracked  under 
Identical  fatigue  load  and  the  super imposition  of  the  profilograms 
of  the  crack  trajectories  (see  footnote  3). 

The  specimen-load  configuration  and  a  sample  crack. 

The  probability  density  of  the  crack  penetration  depth  for  various 

^  =  Tmin^Y*  in  zeroth  approximation.  Parameters  a  and  p  = 

/  -2 
r/£„  are  3  and  10  *  respectively. 

G 

The  probability  density  of  the  crack  penetration  depth  for  various 

A.  Parameters  a.  p  -  r /A  and  q  =  y  4  h*  are  3 .  10~2  and  0. 

u  cm 

respectively  (see  footnote  17). 


Discretization  of  a  crack  trajectory  for  the  purpose  of  the  prob¬ 
ability  calculation. 
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